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MSIBACT 


Face seals In rocket engine fuel and oxldlser turbopumps have been the 
subject of Intense Investigation for over 25 years. Whilst advances have 
been made in the understanding of thin film lubrication between seal faces; 
valuable data has been produced on the friction and wear of material pairs 
In cryogenic environments; pioneering work has been done on the effect of 
lubricant phase change In seals; and many Improvements have been made in 
mechanical seal design, relatively superficial attention has been given to 
the vibrational dynamics of face seals In high-speed turbomachinery. 

BHR Group Ltd. (formerly BHRA) has recently completed the first stage of a 
study, commissioned by the European Space Agency, to Investigate this area. 
This has Involved the development of a two-dimensional adiabatic, turbulent 
lubrication model for thick gas film applications, the production of an 
Integrated mathematical model of gas seal vibrational dynamics for thin 
film applications, Implementation In software, the undertaking of an 
experimental programme to validate software against variations In 
operating conditions and design variables, and suggestions for Improved seal 
design. 
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CURRENT PROJECTS 


• Applications 

rotary 

reciprocating 

static 

• Technology 

experimental 
analytical 
design studies 

• Organisation 
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consortium + Dept, of Trade and Industry 
consortium 4- EEC 
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FLUID SEALING TECHNOLOGY 


Independent facilities and expertise:- 

• Seal Analysis:- 

• Seal testing:- 

• Pump Loops:- 

• Site Measurements:- 

• Design Audit:- 

• Rig Manufacture:- 


thermal, mechanical, lubrication 

oil, water, gas, cryogens, contaminants 

oil, water, slurry, water/air 

fixed, portable 

analysis, critical review 

design, build, modification 
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FLUID FILM 


STRUCTURAL 

DISTORTION 


SEAL RING 
DYNAMICS 



SEAL 

TESTS 


INTEGRATED 

COMPUTER 

MODEL 


Validated by 

J 


PARAMETRIC 

STUDY 

\ 

RECOMMENDATIONS 
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LITERATURE SURVEY 

MODELLING - FILM LUBRICATION ANALYSIS 

Laminar, isothermal - Reynolds equation 

• Liquid films - treatment of cavitation 

• Gas films - grid design (adaptive, graded, etc.) 

- algorithm design (implicit, alternating, 
multigrid, 'interior co-location*, etc.) 

Turbulent lubrication - Hirs* bulk flow 

- Constantinescu 

- Ng & Pan 

Non -isothermal (higher Mach number) compressible flow 

• 1-D (radial) adiabatic model with radial taper 
and entrance effects - Zuk 

Two-phase (boiling interface) films 

• 1-D models 

• Stability approached from consideration of equilibrium 
film thickness vs. opening force curves 

(i.e. not from dynamic analysis) 


LITERATURE SURVEY 
MODELLING - DYNAMICS 

• IK1 and [Cl matricos from fluid film analysis 
then dynamics as a separate problem 

• Integrated analysis - fluid film forces and moments 
calculated at each timestep 

• Excitation mechanisms 

• Number of vibrational degrees -of -freedom 

• Thermal and vibrational transients - 3 or 4 orders of 
magnitude difference in timescale - separate problem 


b H /v,rou P 


bf/A.nn.p 


201 


S©©IPi ©IF MOPilL 

• Concentrate initially on 
DYNAMICS 

• Gas seal assumed 

(2-D transient 2-phase prohibitive 
within commercial constraints) 

• Transient structural distortions 

• Full transient lubrication analysis 

• Turbulent, adiabatic AND laminar, 
isothermal leakage flow 

• Choked exit conditions catered for 

• 4 vibrational degrees of freedom 

• Mechanical damping 
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MATHEMATICAL MODEL 
TURBULENT, ADIABATIC FLOW 


Iterate round 
instantaneous 
equations 


Radial and circumferential velocities: 



2 (optional 


G r and G m from Hirs' bulk flow turbulent 
lubrication theory (or - 1/12, laminar) 

Given in terms of Reynolds numbers "as seen by* 
rotor and stator (different) 

Shear stresses from these Reynolds numbers 

2-D adiabatic energy equation relates pressure 
to circumferential shear stress for density at 
current timestep 


Timestepping 


Substitute u r and u f in continuity equation as 
•Jcnowns' and find density at next timestep 


Use lilm thickness, h, at mid -timestep throughout procedure 
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LAMINAR, ISOTHERMAL FLOW 

Compressible Reynolds equation with ideal gas assumption 


fpfa} ■ 12r t — - {phj 


Time discretisation: 



No energy equation required 

Velocities and shear stresses from pressure gradients 

• Check, exit Mach number distribution for condition 

M < -4= (isothermal flov) 

VT 

• Check Reynolds number distribution for condition 

Re < * e critical (laminar flov) 




)S©IK)AI}OII§(K): 


rotor misalignment 
shaft axial vibration I j 6 


rotor eccentricity 

e 


face coning 


face waviness 









FACE CONING 


Sources: 

• Pro -lapped taper and clamping forces 


Bellows load 



Rotor centrifugal inertia 



• Chang* in ambient temperature for seal ring assemblies 


Sealed pressure differential 


• Interfacial heat generation 

1 

• Interfacial pressure distribution 

*? 












MATHEMATICAL MQDFl.l.fNG 
TIMESTEPPING 


COMPUTATION 

USES 

TO PREDICT 

Fluid film shape 

h *-n i* 


Fluid film pressure 
distribution, p 

p . \.Ij *..'4 

K+x 

Velocity aha shear 
stfets distributions 

K*t 


Coning 

p *‘, 

n*. 

Fluid film forces 
and moments 

M*chabical damper 
forces and moments 

s> *“ '*>.1 \..l 

(X). <*). T * T,or 
Upanxiqi 

in* 


Dynamics 

(XI. (XI. 



Output oY fundamental 
importance 
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MATHEMATICA L MODELLING 
APPROACH TO THERMAL TRANSIENT CONING 

PROBLEM 

~ 

Thermal diffusion timescale 

2-4 orders of magnitude slower than 

timescale of vibrational dynamics 


aMwmow 

1 f 1 1 ""* "' ^ ,l 11 1 '■ 

* Off-line thermal transient Fi. analysis 
provides coning as a function of time 

• Coning-Ume cuve 'sampled* at user-specified 
intervals to provide quasi-steady coning 
for vibrational analysis 
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IT VIPI 


Bellows Damper 



Noseplece Damper 


Leakage Path 


Bellows 


Stator Assembly 


Mating Ring 


Plain-faced balanced mechanical 
seals 

Bellows-type flexibly-mounted 
stator ring 

Rigidly-clamped rotor ring 
Coulomb friction vibration dampers 


bf/Am.up 
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TEST PROGRAMME 

SCOPE 

• 42 tests co vering : 

- face coning 

- rotor eccentricity 

- rotor out-of -squareness 

- degree of damping d 

• Effects of high and low temperature investigated 

• Tests at high speed (60,000 rpm) 

• All tests on typical plain face seal, modified to 
suit required conditions (face diameter « 30mm) 


room temp 
14000 rpm 
5,10.15 bar 


FACILITY 

— 

• Hot gas (to around 220 °C) supply 

• Cold gas supply (polling to room temp) 

• Liquid cryogen supply possible 

• High pressure (rated to 20 bar) up to 14000 rpm 

• High speed (60,000 rpm) at lower pressure 




TEST PROGRAMME 

MISALIGNMENT AND ECCENTRICITY 


MISALIGNME NT 


ECCENTRICITY 



hH /T,nmp 
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• Liquid lubricant film 


• Cavitating film 

• Mechanical contact 

• Circumferential EHD 

• Dffferent spring and secondary 
seal types 

• Floating rotor types 


• Ring seal geometries 





